Stability of de Sitter solutions sourced by dark spinors 
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Cosmology of ELKO and Lorentz Invariant NSS has been investigated using dynamical system 
method starting from the proposal of Boehmer et al [arXiv: 100 3. 3858 . Some different results have 
been obtained by a different approach. The exact solutions described by fixed points of the dynamical 
system and their stability have been discussed. Some stable solutions corresponding to de Sitter 
universes have been obtained. 
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Recently a new class of spinors refer to as dark spinors or ELKO has been found [1] and has received quite some 
attention [2]. Their effects in cosmology have been investigated [3,4,5]. Their dynamics is more general that than 
that of Dirac or Majorana spinors, even when self-interactions are taken into account. This leads to a more general 
and thus more interesting cosmological behavior than that exhibited by normal spinors, including for instance the 
existence of non-trivial de Sitter solutions. However, ELKOs spinors are not Lorentz invariant and their definition 
requires a preferred space-like direction. In a recent article[6], a new class of theories of non-standard spinors (NSS) 
has been constructed. By providing a general treatment of class of NSS models based on an action principle a Lorentz 
invariant, ghost-free, non-local spinor field theory has been developed. The cosmological applications of both the 
original ELKO and the Lorentz invariant NSS have been examined respectively. Especially, and the existence of 
non-trivial de Sitter solutions in each case has been discussed. 

The cosmological dynamics of the effective scalar degree of freedom in both ELKO and Lorentz invariant NSS 
cosmology show a large number of very interesting properties. The cosmological evolution of the NSS energy density 
exhibits a much wider range of behavior than that seen with Dirac spinors. The existence of stable de Sitter solutions 
means that NSS could represent an alternative to scalar field inflation or dark energy. 

In this paper we investigate the evolution of a model of the Friedmann- Robertson- Walker (FRW) universe in both 
ELKO and Lorentz invariant NSS cosmology using the dynamical systems approach [7, 8, 9]. Some results different 
from [6] will be given. 

II. ELKO COSMOLOGY 

Consider a setting in which we only have ELKO spinor fields coupled minimally to gravity, this means we neglect 
all possible interaction terms. We begin by focusing on the background of a flat FRW spacetime with line element: 



X 



ds^ = dt^ 
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A ELKO spinor can be defined by V' = v(0^i with ^ ^ = 1, where ^ is a constant spinor, and then tp can be treated 
as the only dynamical variable cosmological. The energy density p,^, and pressure p^p, are, respectively [6], 
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where H =a /a is the Hubble scalar. The conservation equation P^p +3H{p^ +Pip) =0 implies that the field equation 
for ip is 
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The Friedman equation now reads 
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and the Raychaudhuri equation is 
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The field equations can be formulated as an autonomous system of two differential equations by defining 

u = ttG <^, 
V = VttGcp, 
h = V^H{t), 



and 



and introducing a new time coordinate 



f{v) = (7rG)V(^2), 
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The definition v ~ yirGip here is different from the one given in [6], which leads to different results from [6]. Then, 
the field equations can be rewritten as 
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and 



h = 
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where Ut — du/dr. These equations (5,6) of u and v decouple from h and then define a dynamical flow in a large phase 
volume, the behavior can be analyzed qualitatively by standard techniques from the theory of ordinary differential 
equations. Of particular interest are the exact solutions which define the fixed points of the system and their stability 
around these fixed points. The fixed points of the system can be obtained by solving Ur ~ and Vr ~ for u and v: 



df{v) _ 2vf{v) 



dv (1 - w2) ' 

u = 0, (9) 

which only corresponds to the equation (5.16) in [6]. The equation ^^' = 2f{v)v/{l — v^) can in principle have 
inflnitely many solutions. For example, f{v) = c/ (l — v^j solve this equation for all values of v and in that case we 
would encounter a fixed line. 

The Jacobian of the system (5,6) is 
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where 



dur 2uv r- ju^ + 2f{v) 2^3 
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du (l-i;2) V l-v^ ^{l-v^){u^ + 2f{v)) ' 

dur_ ^ d\f(v) 2 dfjv) u^ + 2f{v) 2u^ + 4f{v) ^ 

dv dv^ (1~V2) dv "" (1-^2) ^ (l-w2)2 ^ 
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At the fixed point (9) we have 
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The Jacobian Ai has the eigenvalues:— 2 a/3 W ^_^.l , 0. For this nonhyperbolic fixed point the linearization theorem 
does not yield any information about the stability of it and therefore, the center manifold theorem is needed. The 
theorem shows that the qualitative behavior in a neighborhood of a nonhyperbolic fixed point p is determined by 
its behavior on the center manifold near p. Since the dimension of the center manifold is generally smaller than the 
dimension of the dynamical system, this greatly simplifies the problem. 
We shift the fixed point to (0, 0) by setting 



V = V + Vc,Vc = —2\/6\ 



1-v^ 
and write the equation (5,6) as, up to second order. 
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has been used. 

The Jacobi decomposition 
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gives the transformation 



u = y, 

V = X •», (13) 
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and then 
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Vr = Ur- (15) 

The system (12) can be written in diagonal form 

Xr = F{x,y), 

Vt = Xx + G{x,y), (16) 



where 
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(x,2/) G K X R, A = ."^ is a negative eigenvalues of the matrix J' and F, G vanish at the origin (0,0) and 

have vanishing derivatives at (0,0). The center manifold theorem asserts that there exists a 1-dimensional invariant 
local center manifold W'^{0) of (16) tangent to the center subspace (the y = space) at 0. Moreover, 1/^^(0) can be 
represented as 

WiO) = {{x, y)eRxR:y = y{x), \x\ < 6} ; y(0) = 0, ?y'(0) = 0, (18) 

for 6 sufficiently small (sec [9], p. 155) where y'{x) — -^. The function y{x) that defines the local center manifold 
satisfies 

y'{x)F (x, y{x)) - \y{x) - G (x, y{x)) = 0, (19) 

for |a;| < S; and the flow on the center manifold VF'^(O) is defined by the equation 

Xr^F{x,y{x)), (20) 

for all x with \x\ < S. 

According to Theorem 3.2.2 in [10], if the origin x = of (20) is stable (resp. unstable) then the origin of (16) is 
also stable (resp. unstable). Therefore, we have to find the local center manifold, i.e., the problem reduces to the 
computation of y{x). 

The condition (18) allows for an approximation of y{x) by a Taylor series at x = 0. Since y{0) — 0, j/'(0) ~ 0, it is 
obvious that y{x) commences with quadratic terms. We substitute 

y{x) = y2x'^ + ysx^ + y^x'^ -\ 

into (18) and set the coefficients of like powers of x equal to zero to find the unknowns 
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means that this fixed point is stable when Wc = or a saddle-node, otherwise. It is important to note that this stability 
is independent of the potential function f{v) — {-kGYV {tp^) . The equation (8) indicate the whole fixed (point) line 
(9) corresponds to de Sitter solutions, however only the point u = 0, w = is stable. 

III. COSMOLOGY OF LORENTZ INVARIANT NSS 

Now we turn to cosmology of the Lorentz invariant NSS spinor tp introduced in [6]. By defining 

and 

the energy density p^ and pressure p^ can be written as 

p^ = vl/ + y($), 

p^ = y'($)$-y($), (21) 

where V'{^) = dV/d^. Then $ and $ can be treated as the dynamical variables cosmological. The field equations 
of them arc 

$ == 2 [$ - y'($)$j , (22) 

$ = -V'($) $, (23) 

The Friedman equation and the Raychaudhuri equation now read, respectively, 

H^='1W + V{^)], (24) 

and 

3ij2 + 2iy=K[y($)-y'($)$]. (25) 

In contrast to [6], we do not search for the de Sitter solutions directly using the condition p^ = —p^. Instead, we 
will discuss the exact solutions of the equations (22) and (23) given by fixed points of a dynamical system and then 
obtain de Sitter solutions from (24). Using 

a 
and (24), (25) the equations (22) and (23) can be rewritten as 



$ = -6 $ \/r [* + Vi*^)] - 3«;F($)$ - -K*$ + ( -K$ - 2 j V'($)$ - ( -K$ - 2 ) *, 
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which have decoupled from H. (25) becomes 
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then we obtain a dynamical system 



with fixed points 
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-K$ + 2J T/($)+ (-K$-2 j F'($)$ = 0, 

e = 0, 
* = -F($), 



-3kV^($) + 3Ky'($)$ - 4y ($) = 0, 

e = 0, 

* = -$v'($), 



3) C: 



The equation (27) indicates that the points B and C correspond to the de Sitter solution 
The stability of these fixed points is described by the Jacobian 
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At B: 






-3«;t/($c)+3«;t/'($c)$c-4F'($c) - 0, 

e, = 0, 



the Jacobian of M has the eigenvalues: -4v/-V''($c), -^-^'(^c) - ^/-^'(^c) - ^"(*c)$c (2 - |k$c), 

-vZ-V^'l^c) + y/-^'($c) - T/"($c)$c (2 - |k$,). 
i) In the case 

V'{<Pc) < 0, (30) 

when 

o<a>e<^,F"(a>e)>o, (31) 

or 

1>c>^,T^"(<I>c)<0, (32) 

the point B is stable, otherwise or $c < 0, it is unstable. However, if 

^"($c) - 0, 



$r = 0, 



the eigenvalues are — 4-^/— y(<I>c), — 2-\/— y ($c), 0, and then B is a saddle-node i.e. it behaves like a saddle or an 
attractor depending on the direction from which the orbit approaches, 
ii) In the case 

V'{<i>c) = 0, 



the eigenvalues are 0, - J-T^"($c)$c (2 - |k$c), J-V"($c)$c (2 - |k$c) 
When 

$c > ^' ^"(*c) > 0, -y"($e)$c [2 - ^K$c ) > 0, 



^ > $, > 0, V"{<Pa) < 0, -y"($e)$c (^2 - ^«:$e ) > 0, 



the point B is unstable, while 



$c > ^' '^'"(^c) < 0, -y"($e)$c (2 - ^K$c ) < 0, 



it is a center. 
When 



B is unstable, while 



it is a center. 
When 



the point B is a center, 
iii) In the case 



when 



B is unstable, while 



^ > $, > 0, F"($,) > 0, -y"($c)$c (2 - ^K<P, ] < 0, 



$c<0,F"($c)>0, 

$c<0,F"($c)<0, 



V'i^c) > 0, 



^r ($,) > y"($c)$c ( 2 - ^K$, ) , 



it is a center. These results indicate that the stability of the fixed point B is dependent on $c, V'i^c) and y"($c) 
but independent of V{^c)- 
At C: 

$c = o,ec = o,*c = o. 



the Jacobian has the eigenvalues: - ^/3kV{0) , - ^/3KV(0) - 2^/-V'{0), -^/3KV{0) + 2^/-V'{0). 
i) In the case 



when 



C is stable. When 



C is stable, while 



it is unstable. When 



y(0) > 0, (33) 



V'{0) > 0, (34) 



V'{0) < 0, -y/3K,V{0) + 2^/-V'{0) < 0, (35) 

V'{0) < 0,-^/3kV{0) + 2v/-V'(0) > 0, 
V'{0) < 0,^^y3KV{0) + 2^/-V'{0) = 0, 



the Jacobin has a zero eigenvalue: —^/3KV(0), — 2-\/3kT^(0), 0, and then C is a saddle-node, 
ii) In the case 

y(o) < 0, 

when 

V'{0) > 0, 
C is a center, while 

V'{0) < 0, 

it is unstable, 
iii) In the case 

1/(0) = 0, 
when 

V'{Q) > 0, 
C is a center, while 

V'{0) < 0, 
C is unstable. The stability of the fixed point C depends on V{0) and ^'(0). 

IV. CONCLUSIONS 

Cosmology of ELKO and Lorentz Invariant NSS has been investigated using dynamical system method. The exact 
solutions described by fixed points of the dynamical system and their stability have been discussed. Some stable 
solutions corresponding to de Sitter universe have been obtained. In the cosmology of ELKO there exists a stable 
fixed point corresponding to a de Sitter universe while in the cosmology of Lorentz Invariant NSS we encounter some 
stable fixed lines corresponding to de Sitter solutions. These results and conclusions, especially the de Sitter solutions 
and their behavior are different from the ones in [6] . 
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